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Abstract 

Charmonium states are expected to be considerably suppressed in the case of 
quark-gluon plasma formation in high-energy heavy-ion collisions. However, a 
robust identification of suppression patterns as signatures of a deconfined QCD 
medium requires a detailed understanding of the "normal nuclear absorption" 
already present in proton-nucleus collisions, where the charmonium production 
cross sections increase less than linearly with the number of target nucleons. We 
analyse the J/ip production cross sections measured in proton-nucleus collisions 
in fixed target experiments, with proton beam energies from 200 to 920 GeV, 
and in d-Au collisions at RHIC, at y/s NN = 200 GeV, in the framework of the 
Glauber formalism, using several sets of parton distributions with and without 
nuclear modifications. The results reveal a significant dependence of the "ab- 
sorption cross section" on the kinematics of the J/?/> and on the collision energy. 
Extrapolating the observed patterns we derive the level of absorption expected 
at Ei a b = 158 GeV, the energy at which the heavy-ion data sets were collected at 
the CERN SPS. 
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1 Introduction 



According to lattice QCD calculations [1], when hadronic matter reaches sufficiently 
high energy densities it should undergo a phase transition to a "plasma" of deconfined 
quarks and gluons (the QGP phase). Considerable efforts are currently being invested 
in the study of high-energy heavy-ion collisions to reveal the existence of this phase 
transition and to study the properties of the new phase in view of improving our 
understanding of confinement, a crucial feature of QCD. The production yields of the 
quarkonium states should be considerably suppressed by "colour screening" if a QCD 
medium with deconfined quarks and gluons is formed in high-energy heavy-ion colli- 
sions [2]. However, already in proton-nucleus collisions the charmonium production 
cross sections increase less than linearly with the number of binary nucleon-nucleon 
collisions. This "normal nuclear absorption" needs to be well understood so that a 
robust baseline reference can be established, with respect to which we can clearly and 
unambiguously identify the signals of "new physics" specific to the high density QCD 
medium. 

So far, the studies of the J/ip suppression patterns determined by the NA50 and 
NA60 experiments at the SPS, from data collected with Pb and In ion beams of 
158 GeV per nucleon, use a baseline reference established on the basis of proton- 
nucleus measurements performed at 400 and 450 GeV, assuming that the energy 
dependencies of the initial and final state "normal nuclear effects", if any, can be 
neglected. It is important to verify if this simple assumption is supported or not by 
existing experimental evidence, looking at the results reported by experiments made 
at different energies, complementing and placing in a broader context the results to 
be obtained by NA60 from proton-nucleus measurements made at 158 GeV. 

In this paper we present a detailed study of some "cold nuclear matter effects" 
affecting charmonium production in proton-nucleus collisions, as well as their energy 
dependence. Section 2 gives a brief historical motivation. Section 3 describes the 
basic framework of our study: how the charmonium production cross sections are 
calculated and how the absorption of the charmonium states in the nuclear matter is 
evaluated. Comparing our calculations to the midrapidity J/V> production cross sec- 
tions measured in proton-nucleus collisions by NA3 (200 GeV), NA50 (400/450 GeV), 
E866 (800 GeV) and HERA-B (920 GeV), presented in Section 4, we observe a signif- 
icant energy dependence of cx^f in this energy range. Simple parametrisations of this 
energy dependence lead to J/ij) normal nuclear absorption rates at the SPS heavy- ion 
energy, i? lab = 158 GeV, significantly larger than at higher energies, as presented and 
discussed in Section 5. 

2 Brief historical motivation 

The NA50 experiment at the CERN-SPS made a detailed study of J/ip and ip' produc- 
tion in fixed-target proton-nucleus collisions with incident protons of 400 and 450 GeV, 
employing six different nuclear targets (Be, Al, Cu, Ag, W and Pb) [3, 4]. Comparing 
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the production cross sections measured at 400 GeV, for instance, to calculations based 
on the Glauber formalism (neglecting nuclear modifications of the parton densities), 
the J/if> "absorption cross section" was determined to be a 3 J^ = 4.6 ± 0.6 mb [3]. 
A similar value, a J J^ = 4.2 ± 0.5 mb, is extracted from a global fit to the 400 and 
450 GeV J /if; / DY cross-section ratios, where Drell-Yan dimuons are used as refer- 
ence. This value has been used by NA50 [5] and NA60 [6] in the studies of the SPS 
heavy- ion data, collected at 158 GeV, assuming that the energy dependencies of the 
initial and final state "normal nuclear effects" , if any, can be neglected. 

An attempt was made [7] to check the reliability of this assumption, exploring 
the centrality dependence of the J /if; / DY cross-section ratio measured by the NA38 
experiment in S-U collisions at 200 GeV per nucleon. This study gave a J J^ — 7.1 ± 
2.8 mb, a value larger than the one derived from the higher energy data. However, this 
result is not satisfactory given the large uncertainty in the value of a J J^ and the need 
to assume that there are no additional nuclear effects between the proton-nucleus and 
the S-U collision systems. We know, in particular, that the if;' state is considerably 
more suppressed in S-U than in p-A collisions [8], and we should not assume that the 
J /ip state is insensitive to the additional mechanisms of if;' absorption. 

Whether or not charmonium absorption depends on the collision energy still re- 
mains an open question. It is long known that E866, at FNAL, observed less J /if; 
absorption at 800 GeV than seen by NA50 in the same xf ~ region. In terms of 
the very simple "a parametrisation" , 

o- p a = cr x A a , (1) 

NA50 reported a = 0.925 ± 0.009 at 400 GeV [3] while E866 obtained values around 
0.95 [9]. 

At the other end of the energy scale, NA3 reported a ~ 0.94 at xp ~ [10]. How- 
ever, it has meanwhile been observed [11] that an absorption pattern generated using 
the Glauber framework, with a certain absorption cross section, leads to significantly 
different a values when the light target is Hydrogen (used by NA3) or Beryllium (used 
by E866). In other words, if NA3 had used Be as the light target, as did NA50 and 
E866, they would have obtained an a value around 0.92. 

A global average of J/if> absorption cross sections was recently reported [12], as- 
suming that measurements collected with different beam particles and energies should 
reflect a single a^ s value despite the observation that some of the values were mutually 
exclusive. 

Knowing the crucial importance of the normal nuclear absorption baseline in the 
interpretation of the J/ip suppression seen in the heavy- ion data and given that the 
charmonium absorption processes may very well depend on the collision energy [13], 
the NA60 experiment collected (in 2004) proton-nucleus collisions at 158 GeV, the 
energy of the Pb and In beams used by NA50 and NA60, with seven different nuclear 
targets (Be, Al, Cu, In, W, Pb and U). In addition, the PHENIX experiment at 
RHIC should soon report measurements of J /if; absorption in y/s NN = 200 GeV d-Au 
collisions with much better accuracy than that available in Ref. [14]. These two sets 
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of results, obtained at very different energies, should help determine the existence of 
an energy dependence and significantly improve our understanding of the mechanisms 
causing the observed nuclear effects in charmonium production. 



3 Basic elements of the calculations 

This section describes the framework we have used to derive the J/ip absorption cross 
sections. The calculations of the charmonium production cross sections are performed 
with the colour evaporation model (CEM) using several sets of parton distribution 
functions (PDFs). It should be emphasised, however, that our results are derived by 
studying how the production yields change from light to heavy targets and, therefore, 
are essentially insensitive to the specific production model and set of PDFs used. We 
performed calculations with "free proton" PDFs and also with PDFs modified by 
the nuclear environment, using several parametrisations of the nuclear modifications. 
The survival probability of the charmonium states traversing the nuclear matter was 
evaluated in the framework of the Glauber model. 

3.1 Charmonium production cross sections 

The charmonium production cross sections used in the studies reported in this paper 
were calculated using the colour evaporation model, described in detail in Ref. [15]. 
In the CEM, the production cross section of each charmonium state, a i: is assumed 
to be a constant fraction, F i: of the total "closed charm" production cross section, 
calculated as the integral over the cc pair mass, to, of the cc cross section, from 
threshold, 2to c , to twice the mass of the lightest D meson, 2tod = 3.74 GeV/c 2 : 



where Xp is the cc longitudinal momentum fraction in the centre-of-mass frame of the 
two colliding hadrons. 

At leading order in perturbative QCD, the cc hadroproduction cross section is 
given by the sum of two partonic contributions, gluon fusion (gg) and quark- ant iquark 
annihilation (qq), convoluted with the parton densities in the colliding hadrons, A and 




2m 



(2) 



B [16]: 




i 



(3) 



\ft(xi,m 2 )f*(x 2 ,m 2 ) + f^xum 2 )f^(x 2l m 2 )]a^m 2 )} 



3 



where X\ and x 2 are the momentum fractions carried by the two partons and y/s NN 
is the centre-of-mass energy of the nucleon-nucleon collision. In our calculations we 
use leading order parton densities, evaluated at scale m 2 = x\ x 2 s NN . 

A basic feature of the CEM is that da cc / dxF dm 2 fully determines the energy and 
momentum dependencies of all the states. The hadronization of the cc pairs into 
charmonium states is nonperturbative, involving the emission of one or more soft 
gluons. Each state requires a different matrix element, condensed in the fractions Fi, 
calculated for each state: J/ip, ip', Xcj, etc. We note that Fj/^, includes both direct 
J/ip production and feed-down J/ip production through Xc tJ radiative decays and ip' 
hadronic decays. 

3.2 Parton densities in the proton and in the nucleus 

Deep inelastic scattering (DIS) and Drell-Yan measurements performed with nuclear 
targets have shown that the distributions of partons in nuclei are significantly mod- 
ified relative to those in free protons. These nuclear modifications depend on the 
fraction of the total hadron momentum carried by the parton, x, on the momen- 
tum scale, Q 2 , and on the mass number of the nucleus, A. While the mechanisms 
governing these modifications are not yet well understood, several groups have pro- 
duced parametrisations, Si(A, x, Q 2 ), that convert the free-proton distributions for 
each parton i, ff(x,Q 2 ), into nuclear ones, ff~(x,Q 2 ), assuming factorisation: 

f l A (x,Q 2 ) = S l (A,x,Q 2 )xf?(x,Q 2 ) . (4) 

Naturally, the nuclear modifications should depend on the spatial location of the 
nucleon inside the target, with the nucleons at the surface being less "shadowed" 
than those in the core of the nucleus [17]. This effect can be ignored in analyses of 
p- nucleus data samples integrated over collision centrality (as done in this work). 

Figure 1 shows several parametrisations of the nuclear modification functions on 
a Pb nucleus for valence quarks (left) and gluons (right), calculated for the scale, 
Q, suitable for quarkonium production calculations with the MRST2001LO [18] or 
CTEQ61L [19] PDF sets. Since the quark and antiquark distribution functions are 
directly probed by the nuclear DIS and Drell-Yan data, their nuclear effects are rela- 
tively well constrained and all parametrisations give similar results. The connection 
between the measurements and the nuclear gluon densities is much more indirect, 
however, relying on the scale dependence of the F 2 structure function and on mo- 
mentum sum rules connecting the momentum distributions of gluons and quarks. See 
Ref. [20] (and references therein) for a recent review of the problem and of the solu- 
tions explored so far. This publication presents the EPS08 model, which uses inclusive 
hadron production data measured at forward rapidity by the BRAHMS experiment at 
RHIC to further constrain the nuclear effects on the gluon densities. Given that those 
measurements might reflect other physical processes (such as gluon saturation in a 
color-glass condensate), it is not at all clear that adding them in the derivation of the 
nuclear PDFs is justified. Therefore, in our study we prefer to place more emphasis 
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Figure 1: Ratio between the valence quark (left) and the gluon (right) distribution 
functions in a nucleon of a Pb nucleus and in a free proton, according to several 
models. 

on the earlier EKS98 [21] parametrisation, the first global analysis of nuclear effects 
on the PDFs. The nDS/nDSg [22] and the HKN07 [23] parametrisations represent 
alternative analyses, also illustrated in Fig. 1, available both at leading order (used 
in our calculations) and at next-to-leading order. 

Charmonium production at fixed-target energies probes x values in the "antishad- 
owing" region, where the parton densities are enhanced in the heavy nucleus. Since cc 
production is dominated by gluon fusion, a good understanding of charmonium pro- 
duction in p-nucleus collisions is presently hampered by the lack of detailed knowledge 
of the nuclear gluon distributions, illustrated by the spread in the curves shown in 
the right panel of Fig. 1. 

Figure 2 shows the influence on the J/ip production cross section per nucleon 
(in the absence of any final state absorption) of using the EKS98 nuclear parton 
distribution functions rather than those of the free proton. The left panel shows that, 
at fixed-target energies, midrapidity charmonium production should be enhanced in p- 
Pb collisions with respect to the linear extrapolation of the pp yields (see Refs. [24, 25] 
for more details). The right panel shows that the "antishadowing" effect turns into a 
"shadowing" effect when charmonium production at forward x-p is considered. 

3.3 Charmonium survival probabilities 

In the framework of the Glauber model, described in detail in Ref. [26], the prob- 
ability that a given charmonium state, generically represented by ip, produced in a 
p-A collision, traverses the target nucleus unbroken by interactions with the nuclear 
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Figure 2: Changes induced by the nuclear modifications of the PDFs on the cc pro- 
duction cross section per nucleon in p-W collisions, using EKS98: as a function of 
y/s NN at x-p = (left) and as a function of x-p at four different energies (right). 



matter, can be calculated as 

SL = J^t = a J d% J J z Pa(6 ' z) S ^ z) ' (5) 

with 

S ihs& z ) = ex P |- dz ' PA(b, z) a^ bs | , (6) 

where b is the impact parameter of the collision (transverse distance between the flight 
path of the incident proton and the centre of the nucleus) and z is the cc production 
point along the beam axis. This "survival probability" depends essentially on the 
nuclear density profiles, pA, and on the charmonium break-up cross section, cr^ bg . In 
our calculations we used Woods-Saxon density profiles with the parameters given in 
Ref. [27]. No charmonium absorption has been considered in Hydrogen nuclei. 

As we will see in detail later on, the absorption cross section crucially depends on 
whether the PDFs are taken to be those of a free nucleon or those of a nucleon in 
a nucleus. For instance, at SPS energies we obtain the same J /if) nuclear absorption 
pattern using proton PDFs and cr^f = 4.5 mb as we do using EKS98 nuclear PDFs 
and a J J^ = 7 mb, as shown in Fig. 3. If the enhancement of the per-nucleon char- 
monium production cross section caused by initial-state antishadowing is ignored, a 
weaker final-state "effective" absorption is sufficient to obtain the same result. 
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Figure 3: Illustration of the interplay between nuclear modifications of the parton 
densities and final-state charmonium absorption. 

4 Measured charmonium production cross sections 

Table 1 lists the experiments that measured the charmonium production cross sections, 
or ratios, used in our analysis. All experiments detected the J/ip through its decay to 
dimuons or dielectrons, with good dilepton mass resolutions and small backgrounds. 



Table 1: Basic features of the experiments providing the charmonium production cross 
sections (or ratios) considered in the present study. 



Experiment 


E lab [GeV] 


Collision systems 


Phase space 


NA3 [10] 


200 


p-H, Pt 


0.0 <x F < 0.7 


NA50 [3] 


400 


p-Be,Al, Cu,Ag,W,Pb 


-0.425 < y cms < 0.575 


NA50 [4] 


450 


p-Be,Al, Cu,Ag,W 


-0.50 < y cms < 0.50 


E866 [9] 


800 


p-Be,W 


-0.10 <x F < 0.93 


HERA-B [28] 


920 


p-C,W 


-0.34 <x F < 0.14 


Experiment 


Js NN [GeV] 


Collision systems 


Phase space 


PHENIX [14] 


200 


pp, d-Au \y cms \ < 0.35, 1.2 < |y cms | < 2.2 



Since our study focuses on the nuclear dependence of the quarkonium production 
cross section, it is preferable to have several target nuclei at each energy. NA50 
collected data with five or six different target materials, from Be to Pb, providing a 
more detailed quarkonium absorption pattern as a function of the size of the target 
nucleus than the experiments which only used two targets. The J /if) production cross 
sections measured by NA50 in p-nucleus collisions, at 400 and 450 GeV, in the centre- 
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Table 2: The J/?/> production cross sections, times branching ratio into dimuons, 
measured by NA50 at 400 and 450 GeV. The errors include statistical and target-to- 
target systematic uncertainties, added in quadrature. 



B x a J/ ^ / A [nb/nucleon] 
NA50-400 NA50-450 "LI" NA50-450 "HI" 



Be 


4.717 ±0.10 


5.27 ±0.23 


5.11 ±0.18 


Al 


4.417 ±0.10 


5.14 ±0.21 


4.88 ±0.23 


Cu 


4.280 ±0.09 


4.97 ±0.22 


4.74 ±0.18 


Ag 


3.994 ±0.09 


4.52 ±0.20 


4.45 ±0.15 


W 


3.791 ±0.08 


4.17 ±0.37 


4.05 ±0.15 


Pb 


3.715 ±0.08 







of-mass rapidity windows —0.425 < y cms < 0.575 and \y cms \ < 0.5, respectively, are 
collected in Table 2. The 450 GeV values correspond to two statistically independent 
data sets, collected with "low" (LI) and "high" (HI) intensity proton beams. The 
400 GeV errors are dominated by a target- dependent relative systematic uncertainty 
of 2.1 %. An extra global systematic uncertainty of 3 %, due to common normalisation 
uncertainties, is not included because it does not affect the evaluation of the nuclear 
dependence. This is because all the 400 GeV p-A data sets were collected in the same 
week, changing the target exposed to the beam roughly every hour, using a rotating 
target holder. On the contrary, the 450 GeV data sets were collected in different 
running periods, over five years, and have independent normalisation uncertainties 
(due to beam counting, trigger efficiencies, etc). Therefore, the quoted 450 GeV 
errors reflect the total uncertainties (added in quadrature). 

The NA50 450 GeV cross sections were also reported in four equidistant xp bins [4] . 
However, given their large uncertainties and small xp coverage, they do not really 
provide extra information with respect to the integrated values, re-analysed in Ref. [3] 
to ensure consistency with the 400 GeV data analysis. Since NA50 reported absolute 
production cross sections (in nb) for each p-nucleus system, we extract the J/ip nuclear 
absorption cross section, cr a j^f, together with a <jq normalisation factor. In the case 
of the other experiments we used a one-parameter fit to extract a J J^ from the ratios 
between the per-nucleon cross sections obtained with "heavy" and "light" nuclear 
targets. The HERA-B and E866 data samples cover a relatively large range in xp 
and are very important for a differential study of the cold nuclear matter effects. 
Table 3 gives the E866 J/-0 W/Be ratios (see Ref. [9] for the xp > 0.2 values) and 
the HERA-B W/C ratios [28] (derived from the exponent a using Eq. 1) which we 
have used in our study. The errors represent statistical and point-to-point systematic 
uncertainties, added in quadrature. The global normalisation errors of 3 % (E866) and 
4% (HERA-B, obtained from a 1.5% uncertainty on a) are not included in Table 3 
but must be considered when comparing different experiments. E866 also measured 
Fe/Be ratios but only for xp > 0.2. 
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Table 3: J/tp cross section ratios measured by E866 and HERA-B, without including 
global errors (3% and 4%, respectively). 



E866 HERA-B 



x-p range 


(x F ) 


W / Be ratio 


x-p range 




(x F ) 


W / C ratio 


-0.10 /-0.05 


-0.0652 


0.8929 + 0.0184 


-0.34 


/-o 


.26 


-0.285 


1.105 + 0.158 


-0.05 / 0.00 


-0.0188 


0.8682 ± 0.0084 


-0.26 


/-o 


.22 


-0.237 


1.034 + 0.096 


0.00 /+0.05 


+0.0269 


0.8720 ± 0.0060 


-0.22 


/-o 


.18 


-0.197 


1.090 + 0.063 


+0.05 /+0.10 


+0.0747 


0.8739 ± 0.0057 


-0.18 


/-o 


.14 


-0.158 


1.043 + 0.042 


+0.10 /+0.15 


+0.1235 


0.8652 + 0.0067 


-0.14 


/-o 


.10 


-0.118 


0.986 + 0.030 


+0.15 /+0.20 


+0.1729 


0.8725 + 0.0100 


-0.10 


/-o 


.06 


-0.079 


0.943 ± 0.022 








-0.06 


/-o 


,02 


-0.040 


0.915 + 0.021 








-0.02 


/+o 


,02 


-0.002 


0.916 + 0.025 








+0.02 


/+o. 


,06 


+0.037 


0.902 + 0.036 








+0.06 


/+o. 


.14 


+0.075 


0.866 + 0.063 



Table 4: J/ip cross section ratios measured by NA3 and PHENIX, without including 
global errors (estimated to be 3 % and 11 %, respectively). 

NA3 PHENIX 



xp range H / Pt ratio y cms range d-Au / pp ratio 



0.0 / 0.1 


1.27 + 0.07 


-2.2 /-1.7 


0.95 + 0.23 


0.1 / 0.2 


1.40 + 0.06 


-1.7 / -1.2 


0.90 + 0.21 


0.2 / 0.3 


1.34 + 0.07 


-0.35 /+0.35 


0.85 + 0.17 


0.3 / 0.4 


1.36 + 0.12 


+1.2 /+1.7 


0.68 + 0.13 


0.4 / 0.5 


1.75 + 0.22 


+1.7 /+2.2 


0.59 + 0.12 


0.5 / 0.6 


2.62 + 0.52 






0.6 / 0.7 


3.58 + 1.81 







Table 4 gives the NA3 H/Pt ratios, for several xp bins, as extracted "by eye" from 
Fig. 2 of Ref. [10]. A common systematic error of 3% [29] is not included. NA3 used 
two targets: protons and platinum. When comparing the NA3 values to the other 
results, we should keep in mind that a proton is not exactly a "nucleus". We also 
analysed the presently available PHENIX d-Au/pp ratios in several rapidity bins [14], 
as collected in Table 4. These ratios have essentially no impact in the results of our 
study, given their large uncertainties. 

It might be worth mentioning that we have not included in our study the mea- 
surements reported by the NA38 and E772 experiments. Later analyses of those data 
sets, made in the framework of the NA50 [7, 30] and E866 [9] experiments, respec- 
tively, revealed that those early results were biased, because of wrongly evaluated 
reconstruction efficiencies (NA38) or acceptances (E772). 
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5 Analysis of charmonium nuclear absorption 



5.1 Cold nuclear matter effects 

Now that we have reviewed the available J/ip production measurements in proton- 
nucleus collisions, we can study them in view of deriving the so-called normal nu- 
clear absorption, described in the framework of the Glauber formalism and quantified 
through the "J/V break-up cross section", cr^f, introduced in Section 3.3. Before 
proceeding, however, we note that this is not the only mechanism affecting the per- 
nucleon production cross sections in p- nucleus collisions. As mentioned in Sections 3.1 
and 3.2, the quarkonium production cross sections crucially depend on the parton 
densities (particularly the gluon densities), which are significantly affected by nuclear 
modifications. These modifications have been taken into account in our calculations, 
employing several parametrisations. Other nuclear effects, such as energy loss, for- 
mation times, etc., are likely to be present and would need to be considered in a 
detailed study of all aspects of quarkonium "cold nuclear effects" . In particular, the 
nuclear charmonium production yields measured by E866 are clearly more suppressed 
at forward xf than at xf = 0, as observed in the W/Be and Fe/Be ratios shown in 
Fig. 4. There is a remarkable change in the suppression pattern at xf ~ 0.25, from 
a relatively flat region around xf ~ 0, where the ip' is more absorbed than the J/V, 
to a forward region where both states show the same strong decrease with increasing 

Xp- 




Figure 4: The xf dependence of the J/ip and ip' production ratios, W/Be (left) and 
Fe/Be (right), per nucleon, measured by E866 at _E lab = 800 GeV (points) and calcu- 
lated assuming only nuclear modifications of the parton densities (curves). 

The curves in Fig. 4 show the expected trends when only considering the nuclear 
modifications on the parton distribution functions, without any final state absorption. 
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It is worth noting the antishadowing enhancement expected in the EKS98 and EPS08 
models at xp ~ 0. 

In the present work, we concentrate on the "midrapidity region" and neglect nu- 
clear effects other than initial-state modifications of the parton densities and final-state 
charmonium absorption, calculated as described in Sections 3.1 and 3.2. As in most 
previous studies of charmonium absorption in nuclear matter, we treat the J/0 as a 
single meson passing through the nuclear medium. However, significant fractions of 
the J/0 yield observed in elementary collisions are due to ip' and Xc decays, 8.1 ± 0.3 % 
and 25 ± 5 %, respectively, as recently determined from charmonium hadropro duct ion 
at x F « [31]. 

5.2 Extraction of cr^J^ from the measurements 

From the equations in Section 3, we can see that the per-nucleon heavy-to-light tar- 
get ratio of the charmonium production cross sections decreases as a function of the 
break-up cross section. This dependence has been calculated respecting the condi- 
tions of each measurement, considering the kinematics window, the collision energy 
and the nuclear matter densities. The calculations were performed using three sets 
of proton parton distribution functions: GRVL0 98 [32]; MRST2001 LO [18]; and 
CTEQ61L [19]. In addition to a control calculation with no nuclear modifications 
of the PDFs, labelled NONE, we employ four models of the nuclear modifications: 
EKS98 [21]; nDSg [22]; nDS [22]; and EPS08 [20]. 



E866 

E lab = 800 GeV 




2 4 



8 10 12 14 
° a bs [ mb l 



Figure 5: Illustration of extracting <7 a b s for J/0 and ■0' production from the E866 
W/Be ratio with and without the EKS98 modifications of the PDFs. 

The rather straightforward extraction of <r a b s for each measurement is illustrated 
in Fig. 5, using the E866 J/0 and ip' W/Be ratios in the range —0.1 < x-p < 0.1. The 
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PDFs effectively cancel in the calculations of the heavy/light ratios, even when nuclear 
effects are considered, since charmonium production is dominated by gg fusion. Given 
that the NA50 data were collected with five or six different targets, we performed a 
global two-parameter fit to all the p-A production cross sections within each data set 
(400, 450-LI and 450-HI). 



Table 5: The J/ip break-up cross sections, cr^f, obtained in the Glauber framework 
described in Section 3. The NA50 400 and 450 GeV values correspond to the rapidity 
windows —0.425 < y cms < 0.575 and \y cms \ < 0.5, respectively. Our study uses 
averages of the two 450 GeV sets. 



Exp. 



a 



NONE 



nDS 



lit [mb] 
nDSg 



EKS98 



EPS08 



NA3 



0.05 
0.15 
0.25 
0.35 



3.77 ±0.98 
5.35 ±0.88 
4.66 ±0.98 
4 Q6 + 151 

^• JU - 1.56 



3.94 ±0.99 

5.46 ±0.88 

4.63 ±0.98 
4 71 +!- 4 9 

^- ' 1 -1.54 



4.27 ± 1.00 
5.85 ±0.89 
5.01 ±0.99 
5 07 + L51 

d - u ' - 1.56 



5.79 ± 1.07 
7.38 ±0.95 
6.18 ± 1.04 

5 81 +L56 

°- OL - 1.61 



7.00 ± 1.12 
8.15 ±0.98 
6.38 ± 1.05 

C- fi r> +1.55 

d - uz - 1.60 



NA50-400 
450-LI 
450-HI 



4.83 ±0.63 
4.51 ± 1.58 
4.82 ± 1.10 



4.74 ±0.62 
4.39 ± 1.58 
4.71 ± 1.09 



4.73 ±0.62 
4.39 ± 1.58 
4.71 ± 1.09 



7.01 ±0.70 
6.89 ± 1.76 
7.17 ± 1.22 



7.98 ± 0.74 
7.93 ± 1.83 
8.21 ± 1.28 



E866 -0.0652 


2.371 


-0.0188 


3.00 1 


±0.0269 


2.90 1 


±0.0747 


2.85 1 


±0.1235 


3.071 


±0.1729 


2.891 


HERA-B -0.158 




-0.118 


0.341 


-0.079 


1.391 


-0.040 


2.111 


-0.002 


2.101 


±0.037 


2.461 


±0.075 


3.521 



0.83 
0.77 
0.73 
0.69 
0.71 
0.67 
0.71 
0.67 
0.72 
0.68 
0.74 
0.70 



-0.83 
-0.77 

2 8S +0 ' 73 
-0.69 

-0.70 

-0.66 

z - ou -0.66 
-0.71 
-0.67 



2.32; 



2.65; 



2.61 



2.31 



+ 0.73 
-0.68 



-0.85 
-0.79 

o «o+ - 75 
°- uz -0.71 
-0.72 
-0.68 

Z - Ud -0.66 
-0.69 
-0.65 
-0.69 
-0.65 



3.01 



3.27; 



2.13; 

1.28 



4.67 
5.39 
4.98 
4.36 
3.95 
3.13 



-0.92 
-0.85 
-0.82 
-0.76 
-0.78 
-0.73 
-0.76 
-0.71 
-0.75 
-0.71 
-0.75 
-0.71 



6.06 

6.20 

5.03 

3.81 

2.98; 

1.91 



-0.98 
-0.90 
-0.85 
-0.79 
-0.78 
-0.73 
-0.74 
-0.70 
-0.72 
-0.68 
-0.71 
-0.67 



1.22 
0.34 
1.18 
1.08 
1.21 
1.10 
1.28 
1.15 
1.51 
1.34 
2.43 
2.02 



0.42 
1.38 
1.99 
1.85 
2.09 
2.96 



-1.22 
-0.42 
-1.18 
- 1.08 
-1.20 
- 1.09 
-1.26 
-1.14 
-1.49 
-1.32 
-2.36 
-1.97 



-1.25 
-0.96 
-1.22 
- 1.11 
-1.24 
- 1.13 
-1.31 
- 1.18 
-1.52 
- 1.35 

z - dz - 1.93 



0.96; 
2.04; 
2.76; 
2.58; 

2.5i : 



0.73 
2.34 
3.68 
4.53 
4.32 
4.28 
4.58 



-1.42 
-0.73 
-1.33 
- 1.20 
-1.32 
- 1.19 
-1.36 
- 1.22 
-1.42 
- 1.27 
-1.65 
- 1.45 
-2.56 
-2.11 



2.23 
3.88 
5.08 
5.46 
4.58 
3.94 
3.59 



-1.52 
- 1.35 
-1.43 
- 1.28 
-1.41 

1.26 
-1.42 
- 1.27 
-1.44 

1.29 
-1.63 
- 1.43 
-2.44 

2.02 



Table 5 collects the cr a ^ s values obtained with the CTEQ61L PDFs, including the 
aforementioned nuclear modifications. Other PDFs give essentially the same results. 
The global uncertainties on the cross-section ratios were propagated into the errors 
on a J J^ using a Monte Carlo procedure; they are common to all the x-p bins at the 
level of the ratios but not at the level of u 3 J^. This way, we can directly compare 
results obtained by different experiments. We can clearly see a correlation between 
the extracted crf/f values and the assumed initial-state nuclear effects on the PDFs 
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(already indicated in Fig. 3): a strong antishadowing effect leads to a large aj^ . 
The most backward HERA-B bins, Xf < —0.14, where a > 1, do not require any 
final-state absorption unless strong antishadowing effects, such as those in EKS98 
and EPS08, enhance the initial-state production rates in heavy targets. Given the 
large uncertainties of the presently available PHENIX ratios, we do not include the 
corresponding a 3 j£ values in Table 5. 



5.3 Dependence of (T 3 J^ on kinematics 

The ratio between the p-W and the p-Be J/ip (and ip') cross sections measured by 
E866, shown in Fig. 4, indicates the existence of additional charmonium absorption 
mechanisms at forward xf, perhaps related to formation time and energy loss effects. 
On the other hand, the measured pattern is remarkably flat in the region xp < 0.25, 
an observation naturally interpreted as meaning that the processes responsible for the 
strong absorption at forward xp axe negligible at midrapidity. This justifies confining 
our present studies to the "midrapidity region" , where the nuclear effects are seemingly 
simpler and "Glauber final-state charmonium absorption", as described in Section 3.3, 
may be sufficient to describe most measurements. Furthermore, the SPS heavy-ion 
data samples have been collected close to midrapidity. Obviously, the flat W/Be 
ratios at "midrapidity" directly translate into an equally flat a(xp), as shown in the 
left panel of Fig. 6. 



0.9 



0.8 



J/y 

HERA-B 920 GeV 
E866 800 GeV 
NA50 450 GeV 



□ EPS08 

• EKS98 

★ NONE 



I 



0.5 



-0.1 



0.1 



0.2 



Figure 6: The nuclear dependence of the charmonium production cross sections mea- 
sured by E866, expressed in terms of a (left) and a 3 J^ (right), without the 3% global 
uncertainty. The left panel also includes the NA50 and HERA-B patterns. 

Our calculations also give rather constant a 3 ^ values in the range —0.1 < xp < 
+0.2, with <7^f ~ 3 mb, when we neglect nuclear modifications of the PDFs, as shown 
numerically in Table 5 and graphically in the right panel of Fig. 6. However, when 
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we use nuclear parton distributions with strong antishadowing, we see that aj^ is 
significantly larger at xp = —0.1 than at +0.2, (see Fig. 6-right). Another observation 
can be made from the xp dependence of the E866 data. If we assume that nuclear 
absorption effects on the J/ip can be effectively described by the Glauber formalism 
with a single a J J^ (ignoring formation times, feed-down contributions, energy loss, 
etc), we see a rather striking increase of a J J^ at forward xp, explicitly shown in Fig. 7 
but already clear in Fig. 4 from the difference between the calculations and the data. 
The observation of a stronger absorption at forward rapidity (or xp) does not depend 
on the nuclear modifications model used. 

The HERA-B measurements also span a relatively large xp range, covering a more 
backward window but having a sizable overlap with E866 around x-p ~ 0. Figure 7 
shows the E866 and HERA-B a J J^ values, obtained using the EKS98 nuclear PDFs, 
as a function of xp and laboratory rapidity, yi a b, calculated from xp accounting for 
the J/0 (which increases with Js NN ). The "midrapidity" (xp < 0.25 or yi ab < 5) 
E866 and HERA-B patterns can be empirically parametrised by the same asymmetric 
Gaussian shape but independent magnitudes. 




Figure 7: The extracted crj^ as a function of xp (left) and yi a b (right), as derived 
from the E866 and HERA-B data. The forward E866 data points (open circles and 
stars) are not included in the fitted curves. The boxes represent the total errors. 

Before continuing, we reiterate our original goal in light of what we have observed 
so far. We seek to determine and quantify possible changes of the "midrapidity J/ip 
break-up cross section" as a function of collision energy. More specifically, we want to 
find out whether it is justified to analyse heavy-ion measurements at 158 GeV using 
the 

^absf va hies derived from proton-nucleus data collected at much higher energies 
(400 and 450 GeV). We have now seen that, in reality, there is no single "midrapidity" 
a tit va lue, unless we neglect nuclear effects on the PDFs. In particular, the E866 data 
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indicate that, for nDSg, EKS98 or EPS08, a 3 J^ drops by a factor of 2 or 3 in the range 
< x-p < 0.25, corresponding to the centre-of-mass rapidity range < y cms < 1. 
Such a non-negligible change in a J J^ with the longitudinal momentum of the J /if) 
indicates that our quest is not an easy one. There are several possible scenarios worth 
considering. For example, it could turn out that all data sets, regardless of collision 
energy, exhibit the same a 3 ^ dependence on y cms , in shape and magnitude. If so, we 
can determine the a 3 J^ suitable for the analysis of the heavy-ion data by integrating 
that universal function in the range < y cms < 1, the NA50 spectrometer window at 
158 GeV. This is likely to result in a different a 3 J^ than that obtained by integrating 
the same function in the 450 GeV window, |y cms | < 0.5. On the other hand, if the 
magnitude of a 3 {^ is energy dependent, the change in the rapidity window could, 
perhaps, partially compensate the change in energy, leading to similar <j 3 J^ values at 
158 GeV in < y cms < 1 and at 450 GeV in \y cms \ < 0.5. It is also possible that the 
shape of cr^(?/ cms ) depends on energy, in which case y cms is not a suitable variable to 
describe charmonium absorption. 

5.4 Dependence of crfj^f on the collision energy 

To evaluate which scenario best describes the measurements requires a global survey 
of all available p-A results, obtained at different energies and in different xp or rapidity 
ranges. 

The a 3 J^ values collected in Table 5 are shown in Fig. 8 as a function of the y cms 
variable, for four different parametrisations of the nuclear effects on the PDFs. The 
E866, HERA-B and NA3 values derived using free proton PDFs (top-left panel) can 
be considered flat in the midrapidity range —0.3 < y cms < 1.0. Therefore, in this 
case (as for nDS) we can evaluate <r^f (y cms =0) by simply fitting each data set to a 
constant. 

The nuclear modifications of the PDFs considerably affect the a J J^(y cms ) pattern 
and a non-trivial y cms dependence is needed to evaluate the cr 3 ^ (y crns =0) values for the 
various data sets. Using EKS98 PDFs, for example, the E866 and HERA-B data sets 
suggest the use of an asymmetric Gaussian function with fi « —0.21 mb, a L ~ 0.37 mb 
and <j r ~ 1.1 mb. Since the NA50 a 3 J^ values correspond to broad rapidity ranges, the 
asymmetric Gaussian shape is weighted by the rapidity distribution of the measured 
J /if) dimuons (those which contributed to the derived <7 a b s value), provided in Ref. [7]. 
However, this convolution has almost no impact on the final result. 

The change in the magnitude of c 3 ^ (y cms =0) with collision energy, y/s NN , can be 
observed in Fig. 9, for free proton PDFs ("NONE") and for the nDSg, EKS98 and 
EPS08 parametrisations of the nuclear PDFs. The corresponding numerical values 

abs 

J/i>, 



are collected in Table 6. To determine the a 3 If relevant for the analysis of the SPS 



heavy-ion results, we extrapolate aj^ (y C ms=Q) down to y/s NN = 17.2 GeV (dotted 
vertical line in Fig. 9), using exponential and linear functions. 

The dependence of a 3 J^ on y cms at y/s NN = 17.2 GeV is shown in Fig. 10, for several 
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Figure 8: a J J^ as a function of y cm s, obtained by considering four alternative 
parametrisations of the nuclear PDFs. The boxes represent the total errors. When 
nuclear PDFs are used, the E866 and HERA-B absorption patterns clearly depend on 



nuclear PDFs, while the values of aj^ (y cms =0) are collected in Table 7. Integrating 
these functions in the NA50 heavy- ion rapidity window at 158 GeV, < y cms < 1, 
weighted by the rapidity distribution of the J/ip dimuons measured in Pb-Pb collisions 
(before acceptance corrections) [33], we obtain the values of a J J^ (0 < y cms < 1) given 
in Table 7. 
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Figure 9: Dependence of cr a ^f (y cms =0) on the nucleon-nucleon centre-of-mass energy. 
The curves represent fits with exponential (solid line with error band) and linear 
(dotted line) functions. 



Table 6: cr^ (y cms =0) values extracted from the five analysed data sets and for the 
nuclear PDFs we have considered, including the free protons case. 

Ex P- vllt (2/cms=0) [mb] 

NONE nDS nDSg EKS98 „ EPSQ8 „, 

NA3 4.71 ±0.66 4.76 ± 0.66 6.78 t L fi 7.82+^ 10.551^ 

NA50-400 4.82 ±0.63 4.74 ± 0.62 5.02 ± 0.67 7.24 ±0.73 8.48 ± 0.79 

NA50-450 4.72 ±0.90 4.61 ± 0.90 4.82 ± 0.95 7.25 ± 1.03 8.36 ±1.08 

E866 2.82 2.53 3.43 1°^ 5.13 lg"g 5.68 ^ 

HERA-B 2.131^ 1.93i^ 2.66lgg 4.35±iS 4.67^ 
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Figure 10: aj^ as a function of !/ cms , extrapolated to 158 GeV considering four 
alternative parametrisations of the nuclear PDFs. 



Table 7: The o" a ^(y cms =0) and crj^(0 < y cuis < 1) values evaluated at 158 GeV 
by extrapolating with an exponential function the values derived from measurements 
made at higher energies, by NA3, NA50, E866 and HERA-B. 



N-PDFs 


^bf(2/cms=°) t mb ] 


^(0 < y cms < 1) [mb] 


NONE 


5.5 ±0.8 


5.5 ±0.8 


nDS 


5.6 ±0.8 


5.6 ±0.8 


nDSg 


7.3 ±0.2 


5.2 ±0.2 


EKS98 


8.7 ±0.7 


7.2 ±0.5 


EPS08 


11.8 ±0.6 


7.5 ±0.4 



5.5 Dependence of crf^f on ^ ne energy 

We now consider a different scenario. If we assume that the J/ip is "broken up" by 
interactions with nucleons while traversing the nuclear target, it is natural to study the 
absorption as a function of the J/^-nucleon centre-of-mass energy, y/s^, a quantity 
that reflects both the nucleon-nucleon centre-of-mass energy, ^s NN , and the energy 
of the J/ip. This energy may be written as a function of the J/ip xp as 



i> n 



\ 



m. 



m. 



(7) 
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where and are the mass and transverse momentum of the J/0- When x-p 
and the dependence is neglected, this expression reduces to 



v /s^ = m M /l + -^ . (8) 
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Figure 11: a J J^ as a function of the J/^-nucleon centre-of-mass energy, determined 
from the fixed-target data sets, with four different nuclear PDFs. The boxes represent 
the total errors. 

Figure 11 shows cr^f as a function of y/s~^ for the fixed-target data sets. The 
centres and widths of the asymmetric Gaussians are exclusively determined by the 
E866 and HERA-B points while the magnitudes are independently fitted for each 
data set. The NA3 values are closer to the others here than as a function of y cms 
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(Fig. 8). However, it is clear that there is also no scaling as a function of y / s^7 v ~ since 

the different data sets result in significantly different magnitudes of cr^f; the E866 
and HERA-B values are too low while NA50 is too high relative to a global average. 
The decrease of a J J^ with increasing beam energy is always seen, regardless of the 
nuclear PDF set, including that of the free proton (NONE). 

The absence of a universal curve indicates that the mechanisms determining the 
J/0 nuclear dependence are not accurately described by the simple "Glauber absorp- 
tion model" we used in our analysis. An improved model is needed to properly explain 
charmonium absorption and its dependencies on collision energy and kinematics, in- 
cluding the extra absorption seen at forward xp. 



6 Summary and outlook 

Previous derivations of the J/0 normal nuclear absorption baseline used to search 
for QGP signals in the SPS heavy-ion measurements, collected at = 158 GeV, 
were based on proton-nucleus data collected at 400-450 GeV, assuming that a 3 J^ is 
a "universal quantity", insensitive to changes related to the collision energy or to 
the rapidity window. This paper presents an improved analysis of charmonium pro- 
duction in proton-nucleus collisions. First, we studied J/0 proton-nucleus data col- 
lected in several fixed-target experiments, covering a broad range of collision energies 
( v /s^7 = 20-40 GeV), as well as d-Au data collected by PHENIX at ^/s^ = 200 GeV 
(presently affected by large uncertainties). Second, we considered nuclear modifica- 
tions of the PDFs, employing several models which consistently indicate initial-state 
gluon enhancement (antishadowing) in the midrapidity region of the fixed-target data. 

We observe that, when the nuclear modifications of the PDFs are taken into ac- 
count, cr^f significantly depends on the rapidity of the J/0, even within a relatively 
narrow midrapidity window. In particular, the J/tp nuclear dependence determined 
by E866 in the window —0.1 < xp < +0.2 only looks independent of x-p if the nu- 
clear effects on the PDFs are neglected. Furthermore, the level of cold nuclear matter 
absorption of midrapidity J/0 significantly decreases with collision energy. While 
the specific numerical values depend on the nuclear PDF sets used, the decrease of 
^abs* with energy is a general feature, independent of any nuclear modifications of the 
PDFs. The observation that depends on the rapidity of the J/0 and on the 
collision energy confirms that the simple Glauber-type absorption model commonly 
used in J/0 suppression studies is insufficient to properly reproduce the available 
measurements. The sought-for "universal quantity" is, after all, a multidimensional 
function without obvious scaling features. 

For the moment, in the absence of a more detailed formalism, we used the data- 
driven cr J J^(y cms ) dependence and a simple extrapolation of the y/s NN dependence to 

evaluate the a J J^ value corresponding to the rapidity window covered by NA50 at 
158 GeV, obtaining the values summarised in Table 7. The a J J^ values obtained in 
this study should be seen as effective ones, incorporating the nuclear matter effects 
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on the directly produced J/ip as well as on the more massive \ c and ip' states. Indeed, 
around 33% of the observed J/ip yields result from decays of Xc and ip' states [31], 
which are expected to suffer stronger nuclear absorption than the directly produced 

Having relaxed the assumption of a universal absorption cross section, we can test 
whether it is justified to study the SPS heavy-ion data, collected at 158 GeV, using 
the 

^abff varue obtained from the higher-energy NA50 measurements. As shown in 
Table 7, if we neglect nuclear modifications on the PDFs we derive er^f (0 < y cms < 
1) = 5.5 ± 0.8 mb, higher than the value used so far in the analyses of the SPS 
heavy-ion measurements (4.2 ± 0.5 mb, also using free proton PDFs). If we use the 
EKS98 parametrisation to model the nuclear modifications of the PDFs, we obtain 
a tbs(® < 2/cms < 1) = 7.2 ± 0.5 mb. It is interesting to notice that this is identical to 
the value obtained by NA50 at 400 GeV in the rapidity window —0.425 < y cms < 0.575 
(7.01 ± 0.70 mb, see Table 5). The drop of <j 3 J^ from y cms = to y cms = 1, in the 
EKS98 case, exactly compensates the increase in a J J^ from E\ &h = 400 to 158 GeV. 
Nevertheless, a quantitative re-evaluation of the level of "QGP melting" in the heavy- 
ion data should be performed, also considering nuclear modifications of the PDFs in 
the beam nucleus (in an x range different from that of the target nucleus). 

Our understanding of the cold nuclear effects on J/ip production, and their energy 
dependence, should significantly improve in the near future, thanks to new NA60 
results based on proton-nucleus data collected at 158 GeV and to new PHENIX results 
based on a large d-Au data set. 
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